Quantum Versus Jahn- Teller Orbital Physics in YVO3 and LaVC>3 
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We argue that the large Jahn- Teller (JT) distortions in YVO3 and LaVC>3 should suppress the 
quantum orbital fluctuation. The unusual magnetic properties can be well explained based on 
LDA+C/ calculations using experimental structures, in terms of the JT orbital. The observed split- 
ting of the spin-wave dispersions for YVO3 in C-type antiferromagnetic state is attributed to the 
inequivalent VO2 layers in the crystal structure, instead of the "orbital Peierls state". Alternative 
stacking of afe-plane exchange couplings produces the c-axis spin-wave splitting, thus the spin system 
is highly three dimensional rather than quasi-one-dimensional. Similar splitting is also predicted for 
LaV03, although it is weak. 

PACS numbers: 71.27.+a, 75.10.-b, 71.70.-d, 75.30.-m 



An essential subject, which is responsible for the rich 
physics in transition-metal oxides, is the orbital de- 
gree of freedom (ODF) and its interplay with the spin, 
charge and lattice degrees of freedom Q, 0, 0- For 
cubic perovskites, quantum orbital fluctuation (QOF) 
would be expected due to the degeneracy in orbital sec- 
tor, and a particular spin and orbital ordered phase 
can be selected by maximizing the energy gain from the 

QOF HHIE0- 0n the other hand > Jahn-Teller (JT) 
distortions will lift the orbital degeneracy, and suppress 
the QOF USE EH. The competing QOF versus JT 
physics has thus been the central issue for the recent dis- 
cussions on tig perovskites. A typical example is LaTi03, 
which shows much reduced ordered moment from that ex- 
pected by mean field theory, and isotropic spin dynam- 
ics W. At one side, the relevance of QOF was empha- 
sized |5j by neglecting the crystal field effects and treating 
the systems essentially as cubic. On the other hand, this 
scenario is questioned by the recent examinations, which 
demonstrated that those unusual magnetic properties can 
be understood from polarized orbital caused by the lat- 
tice distortions 0, [ul 13 E3 ) and that the predicted 
orbital entropy from QOF was not observed pL4| . 

YVO3 and LaV03 are tig perovskites with two lo- 
calized 3d electrons per V. Early experimental results 
showed complicated phases for those compounds 0, 0, 
03. First, LaV0 3 has G-type [H] JT distortion and 
the C-type antiferromagnetic (AF) state below 140K . 
YVO3 does the same between 77K< T <lf 6K, but shows 
the C-type JT distortion, and the G-type AF state below 
77K [16j. Very recently, an unusual magnetic structure 
and dynamic in YV0 3 was reported , especially a split- 
ting of the c-axis spin-wave dispersions was observed for 
the intermediate temperature (77K< T <116K) phase 
of YVO3, where the c-axis lattice dimerization is van- 
ishingly weak 0> 0|. The QOF was thus argued for 
those compounds by neglecting the strong JT distortions 



present (about 2~4%). Theoretically, by treating the sys- 
tem as quasi one-dimensional (ID), the spin-orbital su- 
perexchange model was analyzed, and an "orbital-Peierls 
state" due to the formation of orbital singlet was pro- 
posed 0. However, the spin system can not be regarded 
as ID, while the orbital system is due to the destructive 
interference of the interchain exchange processes 6] . Fur- 
thermore, from this QOF picture it is hard to understand 
the observed large JT distortion and its clear tempera- 
ture dependence It is therefore an interesting and 
challenging problem to judge the underlying physics here, 
especially by first-principles calculations. 

In this letter, we will present firm evidences for the 
crucial role of lattice distortion for YVO3 and LaV03. 
The experimentally observed spin orderings (SO) can be 
systematically explained by our LDA+C/ calculations, in 
terms of the JT orbital. The "unexpected features" of 
spin wave in YVO3 0, namely, 1) splitting of c-axis spin 
wave; 2) \J C \ > \J a b\ {Jc and J a b are c-axis and afe-plane 
exchange couplings, respectively), can be naturally ex- 
plained from the structural point of view. Similar picture 
is also predicted for LaV03. We further point out that 
the theoretical simplification to treat the system as quasi 
ID is lack of firm bases. 

The calculations are done based on the LDA+C/ 
scheme in plane- wave pseudopotential method |2lj . 
For such well defined insulating systems with long range 
ordering, LDA+C/ method typically can give reasonable 
results. The parameter [2^ U c g =3.0eV is used to repro- 
duce the experimental band gaps properly 0. For all 
the discussions, we use the unit cell with a « b « c/V2, 
which includes four V sites (sites 1 and 2 in one layer, 
with sites 3 and 4 on top of 1 and 2, respectively), and 
define the local axes x, y, z as the [110], [110], [001] direc- 
tions of the unit cell. We performed ground state calcu- 
lations for three structures, YVO3 at 65 K, at 100K, and 
LaV0 3 at 1 OK El . Four magnetic structures (i.e., 
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TABLE I: The calculated total energies _E(meV/f.u.), mag- 
netic moment M(^.g/site), and band gap E g (eV), for various 
compounds in different magnetic states. For YVO3 (100K) 
and LaVC>3 (11K), two magnetic moments are given for two 
inequivalent layers. Bold numbers correspond to most stable 
magnetic state. 







FM 


A-AF 


C-AF 


G-AF 


YV03 

(65 K) 


E 


45.9 


27.2 


12.8 


0.0 


M 


1.75 


1.72 


1.70 


1.68 


E a 


0.8 


1.0 


0.8 


1.2 


YVO3 
(100 K) 


E 


16.7 


27.6 


0.0 


19.1 


M 


1.75(1.77) 


1.72(1.75) 


1.70(1.72) 


1.68(1.70) 


E a 


0.8 


1.1 


1.2 


1.0 


LaV0 3 
(10 K) 


E 


38.3 


42.0 


0.0 


23.5 


M 


1.74(1.75) 


1.72(1.72) 


1.70(1.71) 


1.69(1.70) 


E g 


0.7 


1.0 


1.2 


0.9 



FM and A-, C-, G-type AF states 0) were calculated 
for each fixed experimental structures. 

Let us start from the low temperature phase of YVO3, 
which has C-type JT distortion (lif. As shown in Ta- 
ble I, among the four magnetic structures, the G-AF 
state is the most stable state, consistent with experi- 
mental observation and earlier calculations The ob- 
tained magnetic moment (1.68/ib) and band gap (1.2eV) 
are also in excellent a gree ment with experimental ones 
(1.72/Xb and 1.2eV |23jl. demonstrating the validity 
of our approaches. The stabilization of G-AF state can 
be naturally explained in terms of the JT orbital as fol- 
lows. By calculating the occupation numbers, it is easy 
to find that the C-type JT distortion stabilize the C- 
type orbital ordering (OO), where yz, zx,yz, zx orbital 
is occupied for four V sites (1,2,3,4) respectively, and xy 
orbital is occupied for all V sites. In this OO pattern, 
the super-exchange (SE) along the c-axis is expected to 
be AF due to the ferro-orbital chain arrangement of yz 
and zx, according to the Goodenough-Kanamori rules 
(GKR). There are two contributions, however, for the 
SE in afe-plane, the FM one due to the nearly orthogonal 
yz, zx orbital arrangement, and the AF one coming from 
xy orbital. The net coupling will depend on the rela- 
tive polarization of those orbital. The calculated orbital 
occupation numbers {n xy ,n yz and n zx ) for one of the V 
sites (shown in Fig. 1(a)) clearly show that the polariza- 
tion between the yz and the zx states is not so strong, 
and the net magnetic interaction is dominated by the AF 
SE from xy state. As the results, we would expect the AF 
exchange coupling both along c-axis and in the afr-plane. 
This will explain the ground state G-type SO. 

Importantly, the calculated OO patterns and magnetic 
moments are not sensitive to the SO. For the fixed struc- 
ture of YVO3 at 65 K, we always get the C-type OO for 
the four magnetic states. Furthermore, even the occupa- 
tion numbers of each orbital do not change so much for 
different magnetic states, as shown in Fig.l. The situa- 
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FIG. 1: The calculated £29 orbital occupation numbers for dif- 
ferent compounds as function of different magnetic orderings. 
Since the OO patterns are fixed by the structures, only the 
occupations for one of the transition-metal site are shown. 



tion is quite general for all the examples considered here, 
and is also true for LaV03. This clearly suggests the 
crucial role of lattice distortions for those compounds. 
The experimentally observed JT distortions are essen- 
tially important to reproduce the correct magnetic or- 
derings. Giving the facts that all those calculated mag- 
netic states have basically the same OO pattern, and that 
the charge gap E g (see Table I) is much larger than the 
spin excitation energies (typically the order of meV), we 
can approximately decouple the spin degree of freedom, 
and treat it in terms of the Heisenberg model. Then the 
exchange interactions can be estimated by mapping the 
calculated total energies for each magnetic state, E(F), 
E(A), E(C) and E(G), to the Heisenberg model. The 
nearest neighboring exchange coupling constants are then 
given by: 



J c = (1/4S 2 )[£(F) 
J ab = (l/8S 2 )[E(F) 



E(G) 
E(G) 



E(A) 
E(A) 



E(Q] 
E(C)\ 



(1) 



where 5=1 is the moment. For YVO3 at 65K, we got 
J c =7.8meV and J a f>=7.5meV, which is quite isotropic as 
suggested by experments Q , and can be reasonably com- 
pared with the experimental values as shown in Table II. 
Up to this stage, we show that the low temperature phase 
of YVO3 can be well explained by the JT orbital physics. 

For YVO3 at 100K, the experimental G-type JT dis- 
tortion will stabilize the G-type OO, in which the yz and 
zx orbital are occupied alternatively (antiferro-orbital) 
along the c-axis, instead of the ferro-orbital chain be- 
low 77K. For the four magnetic states, we all obtained 
the G-type OO pattern. In such G-type OO state, the 
antiferro-orbital chain along c-axis will favor FM cou- 
pling along the chain, while the exchange coupling in 
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TABLE II: The comparison of calculated and experimental 
spin coupling constants J c and J a b (meV) for various com- 
pounds. The numbers with parentheses are two inequivalent 
values as discussed in the text. The symbols t and J point to 
Ref. |jj and |2j| respectively. 







YV0 3 


YVO3 


LaVOa 






(65 K) 


(100 K) 


(10 K) 


Cal. 


Jc 


7.8 


-7.2 


-6.5 




Jab 


7.5 


0.8(5.3) 


5.8(7.7) 


Exp. 


Jc 


5.7+ 


-2.0(-4.2) T 


-4.0* 




Jab 


5.7+ 


2.6 T 


6.5' 



afe-plane are basically the same as that of YVO3 at 65 
K, resulting in the experimental C-AF state (see Table 

1) [24]]. Now the question is how to understand the "un- 
expected features" of spin wave, which were argued to 
be the result of orbital dimer formation 0. l)The c- 
axis spin wave splits into two branches and open a gap; 

2) |J C | > \J a b\, whil e according to GKR, FM SE (J c in 
this case) is generally weaker than AF SE (J a b)- X-ray 
diffraction results [ljj suggested that YVO3 in the C-AF 
phase has P2\/a symmetry, which has no c-axis dimer- 
ization. The recent far-infrared spectroscopy data 
suggested the possible lowering of the symmetry group 
to P&ll or P\. However, this lowering will not violate 
our follwoing discussions due to: 1) the emergence of new 
phonons is one or two order weaker in intensity |l9l| com- 
pared with the main modes; 2) the following arguements 
are common for all those possible symmetries. A char- 
acteristic point of the lattice structure of YVO3 in C-AF 
state, in sharp contrast with the low temperature phase 
(Pbnm space group), is the absence of any symmetry op- 
eration to transfer one VO2 layer to the neighboring layer 
along c-direction, resulting in two inequivalent VO2 lay- 
ers, which have different amounts of JT distortion. As 
the results, we obtain two J a b (=0.8 and 5.3meV) for 
two different layers [2j|, which stacking along the c-axis 
alternatively. This is in qualitative difference with the 
experimental analysis, in which they assumed alterna- 
tive J c but same J a b to fit the experimental spin-wave. 
By using the calculated exchange parameters (alterna- 
tive J a b), our obtained spin- wave dispersions (shown in 
Fig. 2(B)) definitely shows a c-axis spin-wave splitting, 
which is comparable in size with the experimental one 
(about 5meV). The overall shape of our obtained spin- 
wave is also in good agreement with the experimental 
one. We conclude that the observed spin-wave splitting 
is due to the inequivalent VO2 layers in this compounds. 
For such strongly coupled systems both spin and orbital 
behaviors are essentially three dimensional. Especially 
for the spin degree of freedom, the exchange coupling in 
a6-plane will dramatically affect the spin wave behavior 
along c-axis. It is generally not suitable to treat the sys- 
tem as quasi ID as assumed in previous studies [?]. The 
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FIG. 2: The calculated spin- wave dispersions for various com- 
pounds in the Heisenberg model. Solid lines are results using 
calculated exchange coupling constants, while dotted lines are 
that using single average of two J a b- 



"orbital Peierls states" will be easily suppressed by the 
increased JT distortions. 

The dramatic difference of two J a b comes from differ- 
ent distortions of two layers. For YVO3 at 100K, the JT 
distortion is about 4% for one layer (say layer 1), while 
about 2% for another layer (say layer 2). Such a dif- 
ference will modify the orbital polarization as shown in 
Fig. 1(b) and (c). The polarization between the yz and zx 
orbital for layer 1 is larger than that of layer 2, while the 
n xy is almost same for two layers. The larger polariza- 
tion tend to enhance the FM coupling coming from the 
orthogonal yz and zx occupation in the afr-plane, which 
will compensate the AF coupling from xy orbital. There- 
fore, the net AF coupling are significant suppressed, from 
5.3meV for layer 2 to 0.8meV for layer 1. This fact will 
also explain why | J c | > \J a b\ in this case, and again sug- 
gest the crucial role of lattice distortions. 

Due to the same symmetry, we obtain the similar re- 
sults for LaV03 at 10K, i.e., the observed G-type JT 
distortion will stabilize the G-type OO robustly, and as 
the result, the C-AF ground state is obtained. Our calcu- 
lated magnetic moment (1.7^b) is somehow larger than 
the early experimental value (about 1.3/xb |2oj|). The c- 
axis spin- wave gap is also predicted in this case, although 
it is weak due to the smaller structural difference between 
the two V0 2 layers (see Table I, II and Fig.l, 2). This 
prediction should be verified by future experiments. 

Finally, we have two comments for the C-AF phase of 
YVO3. First, the experimental moment (about 1.05^b) 
is much smaller than our calculated one (1.70 and 1.72 /j,b 
for two inequivalent V sites). This fact has been used as 
an argument for the QOF nature. However, from the 
neutron data , we find that the C-type spin diffraction 
intensity above 77K is far from saturated. By extrapolat- 
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ing that intensity versus temperature data down to OK, 
we can easily get a increased moment by a factor of 1.7. 
Therefore, the reduced magnetic moment could be due 
to the elevated temperature. Second, the spin canting 
of 16 degrees is argued from experimental side Q, while 
it is not included in our calculations. Nevertheless, we 
should say such a canting may quantitatively affect our 
results discussed above (such as the calculated moments 
for YVO3 and LaVC>3 in C-AF state), but not qualita- 
tively, especially for the main conclusion about the struc- 
tural origin of spin wave gap. 

In summary, we show that the reported unusual mag- 
netic properties for YVO3 and LaVC>3 can be reasonably 
explained by a systematic JT picture from LDA+Z7 cal- 
culations. For fixed structures, the obtained 00 patterns 
are not sensitive to SO, suggesting that any meaningful 
orbital fluctuation must be via the phonon degrees of 
freedom. It will be an interesting future subject to study 
the phonon-mediated orbital fluctuations for those com- 
pounds. 
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